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Using 85.2 ±3.6 pb _1 of pp collisions collected at y/s=1.8 TeV with the D0 detector at Fermilab's 
Tevatron Collider, we present the results of a search for direct pair production of scalar top quarks 
(t), the supersymmetric partners of the top quark. We examined events containing two or more 
jets and missing transverse energy, the signature of light scalar top quark decays to charm quarks 
and neutralinos. After selections, we observe 27 events while expecting 31.1 ± 6.4 events from 
known standard model processes. Comparing these results to next-to-leading-order production 
cross sections, we exclude a significant region of t and neutralino phase space. In particular, we 
exclude the t mass m t ~ < 122 GeV/c 2 for a neutralino mass of 45 GeV/c 2 . 
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Supersymmetry (SUSY) 0,0,0, one of the major extensions of the standard model (SM), introduces additional 
particle states. For every bosonic SM particle, it assigns a fermionic "superpartner" and for every SM fermion, a boson. 
The hypothesized SUSY particles include gauginos and scalar quarks or "squarks." The gauginos, superpartners of 
the gauge particles, include neutralinos (prime candidates for dark matter). Squarks include the left-handed and 
right-handed scalar top quarks or top squarks. These weak eigenstates mix to provide the mass eigenstates t\ and ti. 

Generic SUSY searches often make the simplifying assumption of mass-degeneracy of first and second generation 
squarks. The scalar top quark masses, however, are expected to be substantially smaller than those of all other squarks 
El- If sufficiently light, scalar top quarks should be produced strongly at the Tevatron through qq annihilation 
and gluon-gluon fusion with a cross section on the order of that of the top quark 01 ■ According to the next-to- 
leading order (NLO) program PROSPINO [jj, a 100 GeV/c 2 scalar top quark has a production cross section of about 
12 pb, and a 120 GeV/c 2 scalar top quark of approximately 4.2 pb. 

This analysis is sufficiently general that it applies to a broad class of SUSY models. We make no assumptions 
about gaugino unification, but assume that the lightest neutralino Xi is the lightest supersymmetric particle (LSP), 
with conservation of i?-parity guaranteeing its stability. We consider the special case where the scalar top quark is 
light enough that m t - < mi, + mw + m x° an d m f - < nib + m^+, precluding the decays t\ — > bWxi, ti &Xi > and 

ti ~~* (Xi ~^ or Xi* ~* The dominant decay is then t\ — » ex?, yielding an event signature of two jets 

with missing transverse energy {J$t)- We make no attempt to tag the b or c hadrons in jets. 

Characteristics of scalar top quark signal were studied by generating Monte Carlo (MC) events for various com- 
binations of m ( - and m^o , using isajet 0] with its implementation of ISASUSY ^lj • These events were processed 
through a GEANT simulation of the D0 detector, a simulation of the trigger, and the standard D0 reconstruction 
program. 

The major SM backgrounds expected for this signal are multijet events with artificial $ T and vector boson (VB) 
production associated with associated jets. The VB backgrounds include those producing: neutrinos and jets (Z + 2 
jets — ► vv+2 jets and lU+jets, where the W boson decays to a hadronically-decaying r lepton), leptons from VB 
decays that escape detection, or electrons misidentified as jets, pythia [lj| was used to predict the acceptance for 
W '/Z+jet production, while the vecbos 0,0] Monte Carlo generator was used for W/Z + 2 jets events. In each 
case, the calculated cross sections were scaled to match internal D0 reconstruction and acceptance studies for W/Z 
+ n jets. We also used the cross section for ti production measured at D0 0] and the HERWIG generator to calculate 
the acceptance for ti background arising from top quark decays to an undetected charged lepton, a neutrino, and a 
jet. 

The data correspond to an integrated luminosity of 85.2 ± 3.6 pb -1 collected during the 1994 - 1995 Tevatron 
run. The D0 detector consisted of a central tracking system and a uranium/liquid-argon calorimeter surrounded by 
a toroidal muon spectrometer. A detailed description of the D0 detector and data collection system can be found in 
Ref. 0. Events were collected using a trigger requiring two jets, one with Et > 25 GeV and the second with Et 
> 10 GeV, and $ T > 25 GeV, but rejecting events in which the direction of the leading jet and the ~$ T are aligned 
within a polar angle of 14°. Jets are reconstructed offline using an iterative cone algorithm [18j of radius 0.5 in r/ — cf> 
space. A requirement of at least 2 jets with Et>50 GeV, $ T > 40 GeV, and all jets satisfying a difference in azimuth 
Acf>(jet, $ T ) > 30°, guaranteed full trigger efficiency. To suppress VB backgrounds we removed events with electrons 
or muons with £r>10 GeV. 

Multijet backgrounds dominate this sample and arise when mismeasured jets or a misidentified interaction vertex 
induce an apparent ~$ T . Requiring A^(jet, ~$ T ) < 165° eliminates events with jets back-to-back to the $ T . We reduce 
the number of events with poorly measured jet energies by requiring that the Aip between the $ T and the jet with 
the second highest Et exceed 60°. We also removed those events in which jets deposited most of their energy within 
the narrow intercryostat region (0.8 < \i]\ < 1.2) where the central and endcap calorimeters meet. We refer to the 
354 events surviving these criteria as our base sample. (See Table QJ. 

To reduce the background of mismeasured vertices, the central drift chamber (CDC) was used to associate charged 
tracks with jets within the fiducial volume of the CDC, |ry| < 1. Event-by-event these tracks establish the origin of 
each jet, which was required to be no further than 8 cm from the reconstructed event vertex. This vertex confirmation 
was 80% efficient for W — > ev data samples in which electron tracks matched to electromagnetic calorimeter showers 
provided well-defined interaction vertices, while keeping the mis-matched rate below 2%. Table [I] lists the observed 
number of events from the jets plus $ T sample that survive each selection cut down to this clean sample. 

To predict the multijet background remaining in the clean sample, we used events from the base sample where the 
jet vertex position deviated by 15-50 cm from the event vertex. We normalized this background sample to the clean 
sample using events with A0(jet 2, $ T ) < 60° (where jet 2 refers to the jet with the second largest Et). We chose 
the 50 cm value because it provides the best agreement between the background prediction and the data for the f T 
region between 30 and 40 GeV, which is dominated by multijet events. Changing this value to 100 cm (the full width 
of the instrumented interaction region) increases the multijet prediction by 22%, which we take as an estimate of 
the systematic uncertainty of the method. Reversing the order of the vertex confirmation and A0(jet 2,$ T ) selection 
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Selection Events 

2 jets and $ T trigger 536,678 

No detector malfunction or accelerator noise 487,715 

Leading jet E T >50 GeV 205,461 

Second jet £;t>50 GeV 106,505 

# T >40 GeV 13,752 

30° < A0(jet, $ T ) < 165° 4650 

60° < A<t>(jet 2, $ T ) 2327 

Lepton rejection 2009 

All jets reside outside D0 intercryostat region 354 

Vertex confirmation 88 



TABLE I: Number of events surviving in the jets + $ T sample after the application of selection criteria. These 88 events form 
the clean sample. 



Source 


Events in clean sample 


Events in optimized sample 


W/Z 
tt 

QCD Multijet 


63.0 ± 6.9 tu'A 
3.9 ± 0.02 to 2 5 
22.5 ± 7.5 


24.2 ± 3.6 t 9 6 ° 3 
3.4 ± 0.02 t° 2 04 
3.6 ± 1.4 


Total Background 


89.5 ± 14.7 


31.1 ± 6.4 


Data 


88 


27 



TABLE II: A comparison of Standard Model and QCD multijet backgrounds to the number of candidates in the clean and 
final RGS optimized samples. For W/Z/tt the first uncertainty is statistical, the second systematic. For QCD and the total 
observed, the statistical and systematic uncertainties have been added in quadrature. 

with no change in the relative pass rate of events from our base sample showed they provide a legitimate criteria for 
separating subsets for this study. VB background, which include, in decreasing order of importance W — > r + v + 2 
jets, Z — ► v + v + 2 jets and W — > fi + 2 jets, is comparable to the predicted background from multijet production 
(see table |n|). 

A random grid search (RGS) [2(j based on the energy of the two leading jets and the $J T was used to optimize 
the final selection criteria to apply to the clean sample. RGS uses Monte-Carlo-generated scalar top quark events 
to investigate the region of phase space most heavily populated by signal. The RGS was run for the mass points, 
mi = 115 GeV/c 2 and m x ~ = 20 GeV/c 2 , and mi — 130 GeV/c 2 and m Xo = 30 GeV/c 2 optimizing rejection of 
background relative to signal by maximizing the quantity N S i gna \ j J ^Vsignai + -^background ■ This was subject to the 
requirements of > 2% efficiency for signal. While restricting multijet backgrounds to account for no more than 50% 
of the total background. The selection criteria as determined by RGS for each mass point was within 1-2 GeV 
of our final cuts, chosen to be leading jet Et> 60 GeV, second jet Et> 50 GeV, and $ T > 60 GeV. Our final 
sample and estimated background are reported in Table ^ Figure ^ compares data and background for several 
physics distributions. The additional contribution for a 130 GeV/c 2 scalar top quark, 30 GeV/c 2 neutralino signal is 
indicated by the cross-hatched regions on the figure. 

We find the number of observed events is consistent with expected background. Errors on signal efficiencies and the 
fraction of background events passing final selection include the statistical uncertainties from finite MC samples and 
systematics from the jet energy scale (about 7%), luminosity (4.3%) and W/Z cross sections (about 6%). The system- 
atic uncertainty in simulating the trigger is dominated by hardware trigger response (introducing a 5% uncertainty in 
acceptance). The systematic uncertainty in identifying leptons in data ranged from 2-12% (with dependence on the 
detector r\ — <j> position of the lepton and jet multiplicity) and the vertex confirmation procedure includes a systematic 
1% in its efficiency. The signal acceptance of the final selection for a scalar top quark mass of Mj=115 GeV/c 2 , 
M-g =40 GeV/c 2 was 2.7 ± 0.1%. This null result can be represented by a region of exclusion in the (m t ~ , m^o ) plane, 
which is shown in Fig. [21 (along with results from previous experiments). A Bayesian method, using a flat prior for 
the signal cross section and Gaussian priors for background and acceptance, sets the 95% confidence level (CL) upper 
limits. The highest excluded scalar top quark mass value excluded is 122 GeV/c 2 for a neutralino mass of 45 GeV/c 2 . 
The highest excluded neutralino mass excluded is 52 GeV/c 2 for a 117 GeV/c 2 scalar top quark mass. 

We thank the staffs at Fermilab and collaborating institutions, and acknowledge support from the Department of 
Energy and National Science Foundation (USA), Commissariat a L'Energie Atomique and CNRS/Institut National 
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FIG. 1: Data (points) and predicted background (histograms) after final selections. The additional contributions expected 
from a Mj=130 GeV/c 2 , M x ~ =30 GeV/c 2 scalar top quark are shown by shaded histograms. The plots correspond to the Et 
of the leading jet, second jet, $ T (the three parameters optimized using the RGS) and Ht, where Ht =$ T +Y,iET{}eti) , to 
demonstrate agreement with variables not directly optimized via RGS. 
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FIG. 2: Regions in the (mf,m - ) plane excluded at the 95% confidence level assuming 100% branching of t — > cxi- Limits from 

LEP pH and CDF [22^ are also shown in the figure. The dashed lines correspond to kinematic cutoffs from the masses of the 
X?, rn c , mi, and mw- 
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